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In the present paper we study two qubit entanglement in the most general XYZ Heisenberg mag- 
netic chain with (non)homogeneous magnetic fields and the DM anisotropic antisymmetric exchange 
interaction, arising from the spin-orbit coupling . The model includes all known results as particular 
cases, for both antiferromagnetic and ferromagnetic XX, XY, XXX, XX Z, XYZ chains. The con- 
currence of two qubit thermal entanglement and its dependence on anisotropic parameters, external 
magnetic field and temperature are studied in details. We found that in all cases, inclusion of the 
DM interaction, which is responsible for weak ferromagnetism in mainly antiferromagnetic crystals 
and spin arrangement in low symmetry magnets, creates (when it does not exist) or strengthens 
(when it exists) entanglement in XYZ spin chain. This implies existence of a relation between ar- 
rangement of spins and entanglement, in which the DM coupling plays an essential role. It suggests 
also that anisotropic antisymmetric exchange interaction could be an efficient control parameter of 
entanglement in the general XYZ case. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Entanglement property has been discussed at the early 
years of quantum mechanics as specifically quantum me- 
chanical nonlocal correlation Q and it becomes re- 
cently a key point of quantum information theory [J. 
For entangled subsystems, the whole state vector cannot 
be separated into a product of the states of the subsys- 
tems, and the last ones are no longer independent even if 
they are far spatially separated. A measurement on one 
subsystem not only gives information about the other 
subsystem, but also provides possibilities of manipulat- 
ing it. Therefore in quantum computations the entangle- 
ment becomes main tool of information processing, such 
as quantum cryptography, teleportation and etc. 

For realizing quantum logic gates, several models have 
been proposed and demonstrated by experiments in cav- 
ity QED, ion trap, and NMR H, Q. Due to intrinsic 
pairwise character of the entanglement, in all these cases 
important is to find entangled qubit pairs. Generic two 
qubit state is characterized by 6 real degrees of freedom. 
While separable two qubit state has only four degrees of 
freedom. It is clear that single qubit gates are unable 
to generate entanglement in an iV qubit system, because 
starting from separable state we will obtain another sep- 
arable state with transformed by gates separable qubits. 
Then to prepare an entangled state one needs inter qubit 
interactions which is a two qubit gate. The well known 
example of two qubit gate generating entanglement is 
Controlled Not (CNOT) gate 0]. Moreover realization 
of two qubit controlled gates is a necessary requirement 
for implementation of the universal quantum computa- 



tion. For this purpose we also need interacting qubits. A 
simple example of two qubit interaction is described by 
the Ising interaction Jafa2 between spin 1/2 particles. 
More general interaction between two qubits is given by 
the Heisenberg magnetic spin chain model. This model 
have been extensively studied during several decades, ex- 
perimentally in condensed matter systems [1] and theo- 
retically as exactly solvable many body problems (Bethe, 
Baxter and others) [HI- Now they become promis- 
ing to realize quantum computation and information pro- 
cessing by generating entangled qubits and constructing 
quantum gates. Quantum spin chains where proposed 
as medium through which quantum information could 
propogate as a pulse [1^]. Recently several proposals 
discussing quantum entanglement of two qubits in such 
models has been considered [l^. It was noticed that in 
isotropic Heisenberg spin chain XXX model spin states 
are unentangled in the ferromagnetic case (J < 0), while 
for the antiferromagnetic case (J > 0) entanglement oc- 
curs for sufficiently small temperature T < Tc = -j^^- 
Important point is how to increase entanglement in sit- 
uation where it exists already or to create entangle- 
ment in situation when it does not exist. Certainly 
this can be expected from a generalization of bilinear 
spin-spin interaction form. Around 50 years ago to ex- 
plain weak ferromagnetism of antiferromagnetic crystals 
(a — Fe203, MnCOs and CrF^), which has been contro- 
versial problem for a decade, Dzialoshinski (14i | from phe- 
nomenological arguments and Moriya [Tsj from micro- 
scopic grounds have introduced anisotropic antisymmet- 
ric exchange interaction, the Dzialoshinski-Moriya (DM) 
interaction, expressed by 



'Electronic address: 
t Electronic address: 



nilhangurkanQiyte.edu.tr] 



oktaypashaev@lyte.edu.tr 



This interaction arises from extending the Anderson's 
theory of superexchange interaction by including the spin 
orbit coupling effect |15j , and it is important not only for 
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weak ferromagnetism but also for the spin arrangement in 
antiferromagnets of low symmetry. In the present paper 
we show that the Dzialoshinski-Moriya interaction plays 
an essential role for entanglement of two qubits in mag- 
netic spin chain model of most general XYZ form. We 
find that in all cases, inclusion of the DM interaction cre- 
ates (when it does not exist) or strengthens (when it ex- 
ists) entanglement. In particular case of isotropic Heisen- 
berg XXX model discussed above, inclusion of this term 
increases entanglement for antiferromagnetic case and 
even in ferromagnetic case, for sufficiently strong cou- 
pling D > {kTsinh'^e^'^^/^'-'- — J^)^/^, it creates entan- 
glement. These results imply existence of an intimate 
relation between weak ferromagnetism of mainly antifer- 
romagnetic crystals and the spin arrangement in antifer- 
romagnets of low symmetry, with entanglement of spins. 
Moreover it shows that the DM interaction could be an 
efficient control parameter of entanglement in the general 
XYZ model. 



II. XYZ HEISENBERG MODEL 



where 



J X J II 



6^ + J? + and corresponding wave functions are 



1*2) 



1*3) 



1*4) 



V2(/.2 - Bv) 






-{B + y) 

J- 



-{B-y.) 


J++iD 
-{b + v) 




J++iD 




The Hamiltonian of XY Z model for N qubits is 

N-l 

1=1 

{B + b) al + iB-b) af+i + D ■ (a, x a,+i)] 

where B, b- external homogeneous and nonhomogeneous 
magnetic fields respectively, the last term is the DM cou- 
pling. Choosing -j — -j-z the Hamiltonian for two qubits 
becomes 



H = ^[Jx oIgI 4- Jy a\al + a\Gl + {B + b) al 



+ {B-b)Gl^D{alal-alal)\ 
and in the matrix form 



(2) 



For _B = 0, 5 = 0, 13 = the wave functions reduce to 
the Bell states 



*2) - 


|Bo) = 


^ (100) + 111)) 


(3) 


*4) - 


-|Bi) = 


^(|oi) + |io)) 


(4) 


*3) - 


-IS2) = 


^(101) -110)) 


(5) 


*l) - 


-\Bz) = 


^(100) -111)) 


(6) 



The state of the system at thermal equilibrium is deter- 
mined by the density matrix 



H 











2 

J X ~\~ J y 









Jx ~\~ Jy 







+ iD 



J X Jy 

2 




J. 

2 



B 



To study thermal entanglement firstly we need to obtain 
all the eigenvalues and eigenstates of the Hamiltonian ([5]): 
E\^^) = Ei\^i), (i = 1, 2, 3, 4). The eigenvalues (energy 
levels) are: 



El = 
E2 = 



ii 
2 

2 



E3 

Ei 



2 

:k 
2 



PiT) = 



jr^—HjkT-^ 



(7) 



where Z is the partition function, k is Boltzmann's con- 
stant and T is the temperature. Then for Hamiltonian 
^ we find 



-H/kT 



1 + 



-ff 

w 



+ 



2! 



+ ... + 



WkT 



^11 AlA 

A22 A23 

A-i2 A33 

Ail An 



(8) 



3 



where 



^11 

Al4, 

A22 

A23 
A32 
A33 
A41 
A44 



cosh sinh 

kT fi kT 



-e — sinh -f— 
fi kT 



e2kT 



cosh — — sinh 

kT V kT 



^ J++iD . , V 

, 2kT _L smh — 

V kT 

^ J+ - iD . , u 
,2kT -J: smh — 

u kT 



, V b , ^ u 
cosh — — H — smh — — 
kT u kT 



J- 



cosh -— H smh —— 

kT fi kT 



(9) 



and 



Z = Trle'"^'''^] = 2 cosh + cosh -^1 . 

L kT kT J 

As p(r) represents a thermal state, the entanglement in this 
state is called the thermal entanglement. The concurrence 
C (the order parameter of entanglement) is defined as [l^ . 

C = max{Xi - X2 - X3 - X4,0} (10) 
where Xi {i — 1,2,3,4) are the ordered square roots of the 
eigenvalues of the operator 

pi2= pia" <»ay)p*(a^ (»a^) (11) 

and Ai > A2 > A3 > A4 > 0. The concurrence is bounded 
function < C < l.When the concurrence C = 0, states are 
unentangled; when C = 1, states are maximally entangled. 
In our case: 



where p = B , v = \/\P- + D"^ and 



Z = Tr[e 



-H/kT^ 



= 2 



g'jih' cosh + 6^'='-' cosh 

kT kT 



A. Pure Ising Model (B = 0, 6 = 0, D = 0) 

1. Antiferromagnetic Case (J^ > 0): 



The ordered eigenvalues are 



Ai — X2 



> A3 = A4 



-J^/2kT 



where Z — 4: cosh and the concurrence is 

2k 1 



C12 ~ max{ 



JJ2kT 



2 cosh 



-,0} = 



and there is no entanglement. 



2. Ferromagnetic Case (Jz < 0): 



The ordered eigenvalues are 



Ai — A2 



\Jz\/2kT 



> A3 = A4 



\J^\/2kT 



(14) 



(15) 



(16) 



where Z = 4 cosh ttttf and the concurrence is 

_f,-\Jz\/2kT 

Ci2 = max{ il^' 0} — 

2 cosh ^ 

It means that in both antiferromagnetic and ferromagnetic 
cases there is no entanglement in pure Ising Model for any T. 



(17) 



-^3,4 — 



e 2fcT 



Jz 



J- . , n U J- . , U 

1 H 7 smh -7;; =F — smh 

p-' kT p kT 



1 + 



■ sinh -— =F 
kT ^ 



(12) 



smh — — 

kT 



where p = y + Ji , u = ^b'^ + J'^ + D^. Before calculat- 
ing the concurrence for the general XYZ case (Section 7) it 
is instructive to consider particular reductions of XYZ model 
and compare corresponding concurrences with known results. 



III. ISING MODEL 

Let Jx ~ Jy ~ and Jz 7^ and the Hamiltonian is 
= i[Jz + (B + b) al + {B-b) ol + D{olal - ^lol)] 

The eigenvalues are 

A1.2 = ^ (13) 

Jz I 

g2fcT / r)2 ^ J/ D tj 

A3, 4 = — TT- 



, D'^ . ^2 V D . , u 
1 H ^ smh —— =P — smh —— 

1/2 kT^ V kT 



B. Ising Model with Homogeneous Magnetic Field 

(B / 0,6 = 0,D = 0) 



The ordered eigenvalues are 



Ai = A2 



,Jz/2kT 



X3 = X4 — 



-Jz/2kT 



(18) 



where Z = 2 [e^-^-Z^fei' cosh ^ + e''-''^''^] and the 



rence 

C12 = max{ 



,Jz/2kT 



2(e- J./2fcT cosh ^ + eJ-/''T) 
and there is no entanglement. 



.0}=0 (19) 



C. Ising Model with Nonhomogeneous Magnetic 
Field (B = 0,fe 7^ 0,D = 0) 



The concurrence is 
C12 = max{- 



,Jz/2kT 



,0}=0 (20) 



L2(e- J./2fcT ^_ gJ./fcT cosh ^) ' 

and there is no entanglement. 

As we can see in pure Ising model and including homoge- 
neous (B) and nonhomogeneous (b) magnetic fields no entan- 
glement occurs [H, [131, [pl. 
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D. Ising Model with DM Coupling 

(B = 0,& = 0,_D / 0) 



The eigenvalues are 



Ai 



J.h+2D)/2kT 



JJ^-2D)/2kT 



A3 = X4 — 



-Jz/2kT 



(21) 
(22) 



where Z ^2 |^e-^-/2'=T cosh ^ + e'-'-^^ 



1. Antiferromagnetic Case (J^ > 0): 

Ordering the eigenvalues Ai > A2 > A3 = A4 we have the 
concurrence 



sinh 

612 = max{ 



-J,/kT 



cosh ^ + e"-'^^/*'^ 



.0} 



Then C12 = (no entanglement) if sinh < e"''^'""^. When 
sinh > e''^'^'''^ the states are entangled 



C12 



sinh 



-Jz/kT 



cosh + e"-^^/''^ 



(23) 



Moreover states become more entangled for low temper- 
atures: maximally entangled for any D and T — 
so that (limfeT^o C12 = 1) and for stronger DM cou- 
pling limi3_oo C12 = 1. With T growing, Dmin = 
fcTsinh"^ g-^'j/feT jg growing so that we need to increase D 
to have entangled states. 



2. Ferromagnetic Case (J^ < Q): 

a) With weak DM coupling \D\ < |Jz| there is no entangle- 
ment Ordering the eigenvalues A3 = A4 > Ai > A2 we 
have the concurrence 



C12 = max{- 



' cosh 



coshMe-l^.l/fcT + el-^-l/ 



kT 



,0} = (24) 



b) With strong DM coupling \D\ > Ordering the eigen- 
values Ai > A3 = A4 > A2 and the we have the concur- 
rence 



sinh 

(7i2 = max\ 



JJ;.\/2kT 



coshM +elJ=l/2*T 



,0} 



(25) 



IV. XX HEISENBERG MODEL 

For Jz — 0, Jx = Jy = J, the Hamiltonian is 



H=^[J (afaH-a^a|) + (B+6) G't + {B-b) oi+Dialal-alai'^ 



The eigenvalues are 



A1.2 = ^ 



(27) 



A3. 4 = 



1 



1 -I T. smh tt:^ =F smh 



kT 



kT 



where v = \/ J'^ ^ -\- D'^ and 
Z = rr[e-^/'='^] = 2 



,5 , 

cosh — — -I- cosh — — 
kT kT 



A. Pure XX Heisenberg Model (B = 0, 6 = 0, L> = 0) 



The eigenvalues are 

J/fcT 



1 p-J/fcT 

A2 = A3 = -, A4 = ^^. (28) 



1. Antiferromagnetic Case J > 



The ordered eigenvalues are Ai > A2 = A3 > A4 and the 
sinh 7^7 — 1 

concurrence is C12 = max{ , 0} and for 



a) sinh > 1 C12 — 



cosh -j^^T -f 1 ' 

sinh - 
cosh + 1 



so that 



lim C12 = 1 



(29) 



b) sinh if^T < 1 C12 = there is no entanglement for 



T> :i[sinh-il]-i 
k 

Tc 



2. Ferromagnetic Case J < 



(30) 



J^\/kT 



When 



Then C12 = (no entanglement) if sinh -LJ- < e 
sinhM > e^JMkT ^ + *^in(i + e-2|j.|/fcT) ^j^^ 

states are entangled 



sinhM^eU.IAT 
cosh 7^ + el'-'^^l/^^ 

A; J 



(26) 



Moreover states become more entangled for low temper- 
atures limfeT_o C12 = 1 and for stronger DM coupling 
lim_D^oo C12 ~ 1. As we can see there is entanglement even 
in ferromagnetic case with sufficiently strong DM coupling. 
Comparison of (|23|) and (|26p shows that in anti-ferromagnetic 
case, states can be more easily entangled then in the ferro- 
magnetic one 



he eigenvalues are 

^-\J\/kT 

Ai = 



A2 = A3 = 



A4 = 



.\J\/kT 



Z ' ' " Z' 
A4 > A2 = A3 > Ai and the concurrence is 

^ sinh 1^ — 1 
C12 = max{ ,0} 



' cosh ^ + 1 ■ 



and 



a) sinh ^ > 1 



sinh 14 



cosh ^ -I- 1 



(31) 



(32) 
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b) sinh < 1 ^ Ci2 = no entanglement for 



XY HEISENBERG MODEL 



K 



(33) 



In both cases states are entangled at sufficiently small tem- 
perature T <Tc = ^[sm/i"H]"^ 



For Jz = 0, Jx 1^ Jy, the Hamiltonian is 
The eigenvalues are 



B. XX Heisenberg Model with Magnetic Field 

(B / 0, 6 = 0, _D = ) 



The eigenvalues are 

^J/kT 



Ai 



A2 = A3 = — , A4 



■J/kT 



where 



Z = 2 



cosh — — + cosh — — 
kT kT 



and the concurrence is C12 = max{ 



sinh -jyp — 1 
cosh -pj: + cosh ' 



(34) 

(35) 
0} and 



a) sinh ^ > 1 C12 



sinh ^ - 1 
cosh -pp + cosh ^ 



b) sinh -jj^t < 1, C12 = no entanglement for 



. , o u J- . , n 
1 H smh — =F — smh — 

/i2 kT fj. kT 



(39) 



sinh ■ 



kT 



where = ^ B'^ + Jl, 1^ = + 1)2 + J2 , J± = iiiii 
and 

^ = rr[e-^/''^] = 2 [cosh ^ + cosh ^] . (40) 



A. Pure XY Heisenberg Model (B = 0, & = 0, D = 0) 



The eigenvalues are 



T>M[sinh-M]-^ 



(36) 



It shows that inclusion of magnetic field does not change 
the critical temperature for concurrence in both anti- 
ferromagnetic and ferromagnetic cases. Pairwise entangle- 
ment in A'^- qubit XX chain and experimental realization of 
XX model has been discussed in d, [|, [H], [13 and [H. 



Ai = 



where 



A2 = 



Z = 2 



-J^/kT 



A3 = 



A4 = 



-J+/kT 



1 J~ T J+ 

cosh h cosh — 

kT kT 



(41) 



For Jx = J(l + 7) and Jy — J(l — 7) so that J+ — J, J_ 
J7, the eigenvalues are 



C. XX Heisenberg Model with DM Coupling 

(B = b = 0, I> / 0) 



Ai = 



oJl/kT 



A2 = 



-J~l/kT 



A3 = 



,.J/kT 



A4 = 



-J/kT 



(42) 



The eigenvalues are 



1 g-/3/fcr 



A4 



pfi/kT 



(37) 



where /? > 0, /3 = V + 1)2 ^nd ^ = 2(1 + cosh ^). The 
ordered eigenvalues axe A4 > A3 > Ai = A2 and the concur- 
rence is 



, sinh 7^ — 1 
C12 = max{— -i^— ,0} 
cosh 1^ + 1 



(38) 



where !^ = \/7P + TP : 

sinh 7^ — 1 

a) sinh^>l^a2^ ^^^^|^^ 

b) sinh -jT^T < 1 C12 = there is no entanglement. 

Entanglement increases with growth of DM coupling in both 
anti-ferromagnetic and ferromagnetic cases. 



1. Anti-ferromagnetic Case Jx > and Jy > 

The ordered eigenvalues are A3 > Ai > A2 > A4 and the 
concurrence 

„ ^ sinh ^ — cosh , , , 

C12 = max{ f- , 0} (43) 

cosh + cosh -j^ 

J , T sinh — cosh 

a) sinh^ > cosh^ C12 - ''^ ''^ 

(limT^o C12 = 1) 

b) sinh < cosh => C12 = there is no entanglement. 

In Fig. 1, we plot the concurrence C12 in XY Heisenberg 
antiferromagnet as function of ^ and ^ 



cosh + cosh ^ 
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FIG. 1: Concurrence C12 in XY antiferromagnet as function 



of ¥t and ^ 



B. XY Heisenberg Model with DM Coupling 

(B = 0,6 = 0,D / 0) 



The eigenvalues are 



where 



Ai 
A3 



Ao = 



-J^/kT 



-^Jl\+D^lkT 



Z = 2 



cosh ——T + cosh ■ , „ 
kT kT 



(47) 
(48) 

(49) 



1. Antiferromagnetic Case 



2. Ferromagnetic Case Jx < and Jy < 



C12 = 



a) sinh > cosh ^ 



sinh ^^tJ- — cosh 

,0} 



kT 



cosh ^-jpjt + cosh 



sinh ^ — cosh ^-4tJ- 

kl kl 
<-^12 = 7J—, J— 

cosh ^^TjT^ + cosh ^ 



lim C12 = 1 



(44) 



(45) 



(46) 



The concurrence is 
- max{ 



, sinh ^ t„ — cosh 
C12 = maa;| , ''^ ^ 

cosh -^^ — jjTp h cosh -pp 



.0}. (50) 



coHU - — jjTp h cosn jjTjT 

It shows that for any temperature T we can adjust sufficiently 
strong DM coupling D to have entanglement. 



strong, 

a) sinh ^ > cosh 



sinh ^ — cosh 

r~i kl ^ ' 



cosh 



(51) 



kT 



+ cosh ^ 



b) sinh - — jyp — < cosh C12 ~ there is no entangle- 

ment. 



b) sinh ^^7=^ < cosh ^ 



C12 = there is no entanglement. 



In Fig. 2, we plot the concurrence C12 in XY Heisenberg 
ferromagnet as function of ^ and It^- 




FIG. 2: Concurrence C12 in XY ferromagnet as function of 



Tt and l^-l 



Thermal entanglement in XY chain was studied in [2J|, 
[2^ and [25I ] in the presence of external magnetic field B and 
in [231 by introducing non-uniform magnetic field b. 



2. Ferromagnetic Case 

Ferromagnetic case gives the same result as anti- 
ferromagnetic case. Comparison with pure XY model (|43p 
and (|44p shows that level of entanglement is increasing with 
growing DM coupling D and C12 = 1 when D 00. 



VI. XXX HEISENBERG MODEL 

For Jx = Jy = Jz = J, the Hamiltonian is 



H = ^[J{a^ia^2+'TM + '^t'^i) + iB + b)at 
+ {B^b)a'2+D{a^,ay-afa^)]. 
The eigenvalues are 



(52) 



Al,2 = 

A3, 4 = 

where 
Z ^2 



-J/2kT 



^J/2kT 



(53) 



H- 



J-^+D'^ 



sinh — =F 
kT 



sinh ■ 



kT 



-j/2kt , B j/2fcT , yj^TPT^ 

' cosh h e cosh 



kT 



kT 



(54) 
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A. Pure XXX Model (B = 0, 6 = 0, _D = 0) 



2. Ferromagnetic case (J < 0): 



The eigenvalues are 

-^1,2 = 7^ — , As = 

where 



-,//2feT 



z 

Z = 2 



A4 = 



3J/2feT 



-J/2fcT , J/2feT V 

e + e ' cosh -r— 



(55) 
(56) 



Ci2 = max{- 



' cosh 



I , , 0} = and no en- 

coshJfii + el-'IAi'cosh^ 

tanglement occurs. Therefore inclusion of magnetic field does 

not change the result. Entanglement in XXX Heisenberg 

model with magnetic field has been studied in 



Antiferromagnetic case (J > 0): 



C. XXX Heisenberg Model with DM Coupling 

(B = 0, 6 = 0, L» 7^ 0) 



The concurrence is 



g2.7/feT _ 3 

Ci2=max{-——,0} 



a) sinh^ > e-^/'^^ 



Cl2 — 



sinh - 



-J/kT 



e-J/2*T + cosh: 



(57) 



(58) 



2kT 

For sufficiently small temperature T < entangle- 
ment occurs and limT^o C12 = 1 
b) sinh < e~''^^'^ C12 = there is no entanglement. 



The eigenvalues are 

g~J/2fcT 
Al,2 ~ 7^ , 



A3 = 



where 



A4 = 



Z = 2 



JJ+2^ .J^+D^)/2kT 



-JI2kT , JI2kT , \/j2~f~D2 

' + e ' cosh — 



(62) 
(63) 

(64) 



Ferromagnetic case (J <0): 



— cosh 

The concurrence C12 ~ max{ r-^^ — , 0} = 

cosh + el-^l/''^ 

and no entanglement occurs. 

Thus, for ferromagnets, spins are always disentan gled , 
while entanglement is observed for antiferromagnets [9|. [13||. 



B. XXX Heisenberg Model with Magnetic Field 



The eigenvalues are 



Al,2 = 



-J/2kT 



A3 = 



-J/2kT 



A4 = 



3J/2feT 



Z = 2 



~ J/2feT 



^ B 

cosh h e 

kT 



J/2kT 



cosh ■ 



kT 



(59) 
(60) 



1. Antiferromagnetic case (J > 0): 



The concurrence 



C12 ~ max{ 



„2J/kT 



-3 



e2J/fcT^_l4.2coshT|T 



,0} (61) 



a) sinh^ > e--'/*^, C12 



g- j/2fcT cosh ^ + cosh 2^ 



For sufficiently small temperature T < -j^^ entangle- 
ment occurs and limT^o C12 ~ 1. Comparison of (|61[) 
with (|57p shows that inclusion of the magnetic field B 
does not change the critical value but decreases the level 
of entanglement, 
b) sinh < e~''^'^'^ C12 = there is no entanglement. 



1. Antiferromagnetic Case (J > 0^; 
The concurrence is 



C12 = max{- 



, sinh 



kT , V£!±£! 



,0} (65) 



g-j/kT _|_ eosh 



a) sinh 



C12 



sinhV^_g-^A^ 



g J/fcT _|_ gQgJj 



For a given temperature, when 



D > VkTsinh-'^e-J/'^T _ J2 
there is entanglement. 



(66) 



(67) 



b) sinh < e-^/*^ ^ C12 = there is no entangle- 

ment. 



2. Ferromagnetic Case (J <0): 
The concurrence is 



C12 = max{ 



sinh 



„^iji/fe 



jt J|/fcT cosh 



(68) 



a) sinh 



Vj^'+d' > gUI/ 



sinh^-el^l/- 
O12 = 



jlJI/fcT _|_ 



cosh ^ ^ 



(69) 
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For a given temperature, when 

D > \/kTsinh-^e\J\/''T _ ,p 
there is entanglement. 



(70) 



b) sinhv5£!<eUIAT^^^^^Q 

As we can see inclusion of DM coupling D in XXX case 
increases entanglement in antiferromagnetic case and even 
create entanglement in ferromagnetic case. Thermal enat- 
nglement and entanglement teleportation in XXX Heisen- 
berg chain with DM interaction has been studied in 



b) sinh -gp <e ■'"^^'^ ^ C12 = 0, no entanglement. 

From above formulas follow that for sufficiently small 
temperature T the states are entangled. 

For A < — 1 the concurrence is 



C12 = max{ 



- cosh 



coshi^ + elJ-l/'=i' 



,0} =0 



(76) 



and no entanglement. 

For A = 1 the anisotropic model reduces the isotropic 
XXX model, and the concurrence reduces to 



VII. XXZ HEISENBERG MODEL 



C12 = max{ 



„2J/fcT ^ 



g2J/fcT ^ 3 



0} = 



(77) 



For Jx = Jy ~ J Jz the Hamiltonian is 



H = ^[J{ata^ + alay,+Aata^) + {B + b)at (71) 



+ (B - 6) cr| + D(o-Ja| - a^a|) 
where A = Jz/ J- The eigenvalues are 



Al,2 

-^3,4 



e 2fcT 
Z 
e2kT 



(72) 



1 H smh — — =F smh 



kT 



kT 



where fj, = B, v = \'J^+D^ + W and 
Z = 2 



J^/2'=^cosh-^+e^^/^''-^cosh^ 



kT 



kT 



When the temperature is larger than the critical tem- 



perature Tc 
pears. 



the thermal entanglement disap- 



2. Ferromagnetic Case (J < 0): 



For A < 1 The concurrence is 

sinhHl-e'^-l/'=^ , 

C12 = 7nax{ tfj- ,0} 

coshiil +el-^-l/'=T 



a) sinh HI > el-'- 1/'=^ 



sinhm_el-^-IA" 
coshH, +elJ-l/'=T 



(78) 



(79) 



A. Pure XXZ Heisenberg Model 

(B = 0,& = 0,D = 0) 



In this case the eigenvalues become 



Al,2 = 



- Jj/2fcT 



A3 = 



„(Jj-2J)/2fcT 



A4 = 



„(.7^+2J)/2feT 



where /3 = J and Z = 2 ^e'-''/^'''^ + e'''^^'''^ cosh . 



(73) 



1. Antiferromagnetic Case (J > 0): 

For |A| < 1 weak anisotropy (A > easy axis, A < easy 
plane) and A > 1 strong anisotropy, the concurrence is 



sinh ^-6-^-/'=^ , 
C12 = maa;{— ;--^^_____,o} 



'coshH, + e-././fcT' 



a) sinh > e ^ 



J^/kT 



C12 



J 

kT 



-J^/kT 



cosh H, _|_ f^-JJkT 



(74) 



(75) 



b) sinhHl <el-^-l/'=^ 



■ C12 = 
For A > 1 the concurrence is 



— cosh 

C12 = max{ ,0} = 

Kl 



and no entanglement. 



(80) 



B. XXZ Heisenberg Model with DM Coupling 

(B = 0,6 = 0,D / 0) 



The eigenvalues are 

„-J,/2fcT 



A3 = 



„(./,-2^72+I32)/2fcT 



where 



A4 



Z = 2 



e ' + e ' cosh 



+ L>2 



kT 



(81) 
(82) 

(83) 
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1. Antiferromagnetic case (J > G): 
The concurrence is 



Ci2 = max{- 



. sinh 



cosh 



-,0} 



a) sinh 



sinhA^4±^-e-^^/^^ 

kT 



cosh 



(84) 



(85) 




FIG. 3: Concurrence in XY Z antiferromagnet as function of 
T 



b) sinh 



Ci2 = Comparison with ffSl 



shows that with growth of D entanglement increases. 



2. Ferromagnetic Case (J < Q): 



This formula shows that entanglement increases with low- 
ering temperature. In Fig. 3, we plot function /(T) for 
{Jz,Jy,Jx) = (3,2,1). It shows entanglement for T < 
Tc- In addition, from (|90p we have entanglement increas- 
ing with growing anisotropy J+ and decreasing with growing 
anisotropy J_. Moreover it increases with growing Jz. 



a) For small D < Dc — \/ J'^ — J'^ no entanglement. 

b) For D > Dc the concurrence is 



2. Ferromagnetic Case 



Ci2 = max{ 



,sinhVZ!p!_eU.IAT 



cosh^^^+el-'-l/fe^ 



. 0} (86) 



and entanglement increases with growing D. 

Entanglement for XX Z Heisenberg model was considered in 
fH and effect of DM interaction on XX Z model in [H. 



VIII. XYZ HEISENBERG MODEL 



A. Pure XYZ Model (B = 0, 6 = 0, D = 0) 



The eigenvalues are 

(-J,-2J_)/2fcT 

Ai = ^ , 



„(J.-2J+)/2fcT 



A2 = 



„(-J^+2J_)/2feT 



„(Jj+2J+)/2fcT 



where 



Z = 2 



'^^/^''■^ cosh 4 cosh 4 



kT 



1. Antiferromagnetic Case 



kT 



(87) 



(89) 



Let Jz < Jy < Jx < then Ja 



U+l, J- 



U-l < 



and = — I Jz|. The biggest eigenvalue is 

(\J^\ + 2\J_\)/2kT 

Ai = ^ 



and the concurrence is 



^ sinh ^-TTfr 
Ci2 = max{ 



cosh ^ 



l-'+lg-l^^l/fcT 



cosh 



J- 1 



+ cosh^-^e-\J-\/'<^T 



,0}. 



Then entanglement occurs when 



f(T) = sinh 1^ - cosh l^e" 
■'^ ' kT kT 



\J^\/kT 



> 0. 



(91) 



This formula shows that entanglement increase with low- 
ering temperature. In Fig. 4, we plot function f{T) for 




Jz > Jy > Jx > ^ J+ > 0, J= - \J- \ < 0. The biggest 

|./^I+2|.J+I 

eigenvalue is A4 = — and the concurrence is 



C12 = max{ 



sinh ^ — cosh •''/'^'^ 

kl kl 

cosh ^ + cosh x^e"''^'"''^ 



0}. 



Then entanglement occurs when 



/(T) = sinh ^ - cosh ^e-""' > 0. 



-Jz/kT 



(90) 



FIG. 4: Concurrence in XYZ ferromagnet as function of T 

{Jz, Jy, Jx) = (—3, —2,-1). It shows that entanglement in- 
creases with growing anisotropy J+ and decreases with grow- 
ing anisotropy J_. Moreover it increases with growing Jz- 
Thermal entanglement in pure XYZ model has been studied 
in [3J, [2^. Enhancement of entanglement in XYZ model in 
the presence of an external magnetic field considered in [33 |. 
and influence of intrinsic decoherence on quantum teleporta- 
tion in [3311. 
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B. XYZ Model with Magnetic Field 

(B = 0,& = 0,_D = 0) 

The full anisotropic XY Z Heisenberg spin two- qubit sys- 
tem in which a magnetic field is applied along the z-axis, was 
studied by Zhou et al. The enhancement of the entanglement 
for particular fixed magnetic field by increasing the z- compo- 
nent of the coupling coefficient between the neighboring spins, 
was their main finding. 



XYZ Model with DM Couphng 

(B = 0,fe = 0,_D / 0) 



2. Ferromagnetic Case : 
The concurrence is 

sinh 7^ — e' -^'l'"'^ cosh 

Ci2 = max{ ^ , 0} 

cosh -pj; + e^J^l/i^T (,Qgjj __ 

and entanglement occurs for sufficiently strong D 

J- 



sinh 



> e 



cosh ■ 



kT kT' 
Fig. 6, shows Ci2 as function of D and T. 



(96) 



Al 



As = 



„(-.J^+2J_)/2feT 

z ' 

„(J^+2i/)/2fcT 



A2 



,(-J^-2,7_)/2feT 



(92) 



Z 

A4 = ^ (93) 




where v ■ 



Z = 2 



-^^''^'^^ cosh ^+6'^^/^'=^ cosh ^ 

kT 



V 



(94) 



FIG. 6: Concurrence C12 in XY Z ferromagnet as function of 
D and T 



1. Antiferromagnetic Case 
The concurrence is 



sinh-T^-e ■^^^'"'^ cosh ^ , , , 
Ci2 = max{ 1^,0} (95) 

cosh -pp + g- Jj/fcT gQgJj __ 



and entanglement occurs when 



sinh.4^>e-/-cosh§. 

In Fig. 5, we plot the concurrence C12 as function of D 
and T. Comparing with pure XYZ case (|90p . we find that 
inclusion of DM coupling increases entanglement. 
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FIG. 5: Concurrence C12 in XYZ antiferromagnet as function 
of D and T 
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IX. CONCLUSION 

We found in general if Ai (the largest eigenvalue) is degen- 
erate with A2 then no entanglement occurs. From our con- 
sideration follows that in all cases decreasing of temperature 
increases entanglement, if it exists . So that at zero temper- 
ature T = states are completely entangled C12 = 1. This 
fact links entanglement with the Mattis- Lieb [l^ theorem on 
absence of phase transitions in one dimension at T 7^ 0. More- 
over, inclusion of the DM coupling always increases entangle- 
ment, this is why it could be an efficient control parameter of 
the entanglement. Our results show existence of intrinsic re- 
lation between weak ferromagnetism of mainly antiferromag- 
netic crystals and spin arrangement in (anti)ferromagnets of 
low symmetry with entanglement. 

Very recently thermal entanglement of a two-qubit isotropic 
Heisenberg chain in presence of the Dzyaloshinski-Moriya 
anisotropic antisymmetric interaction and entanglement tele- 
portation, when using two independent Heisenberg XXX 
chains as quantum channel, have been investigated [S^. It 



was found that the DM interaction can excite the entangle- 
ment and teleportation fidelity . As was noticed DM inter- 
action could be significant in designing spin-based quantum 
computers [34l ]. Moreover, studying the effect of a phase shift 
on amount transferable two-spin entanglement in a spin chain 
[ssj . it was shown that maximum attainable entanglement en- 
hanced by DM interaction. 

Therefore would be interesting to consider most general 
XYZ Heisenberg models with DM interaction as quantum 
channel for quantum teleportation which requires to know de- 
pendence of pairwise entanglement on the number of qubits in 
the spin chain. These questions now are under investigation. 
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